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In general relativity, all black holes in vacuum are described by the Kerr metric, which has 
only two independent parameters: the mass and the spin. The unique dependence on these two 
parameters is known as the "no-hair" theorem. This theorem may be tested observationally by 
using electromagnetic or gravitational-wave observations to map the spacetime around a candidate 
black hole and measure potential deviations from the Kerr metric. Several parametric frameworks 
have been constructed for tests of the no-hair theorem. Due to the uniqueness of the Kerr metric, 
any such parametric framework must violate at least one of the assumptions of the no-hair theorem. 
This can lead to pathologies in the spacetime, such as closed timelike curves or singularities, which 
may hamper using the metric in the strong-field regime. In this paper, I analyze in detail several 
parametric frameworks and show explicitly the manner in which they differ from the Kerr metric. 
I calculate the coordinate locations of event horizons in these metrics, if any exist, using methods 
adapted from the numerical relativity literature. I identify the regions where each parametric 
deviation is unphysical as well as the range of coordinates and parameters for which each spacetime 
remains a regular extension of the Kerr metric and is, therefore, suitable for observational tests of 
the no-hair theorem. 
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I. INTRODUCTION 

According to the no-hair theorem, the exterior gravita- 
tional fields of isolated and stationary black holes in gen- 
eral relativity are uniquely characterized by their mass 
M and (the magnitude of their) spin angular momen- 
tum J . In particular, this field is described by the Kerr 
metric, the unique stationary, axisymmetric, asymptoti- 
cally flat, vacuum solution to the Einstein equations that 
possesses an event horizon and is free of closed timelike 
curves outside of the horizon [iHal ■ The no-hair theorem, 
therefore, implies that all multipole moments of the Kerr 
spacetime are completely determined by the first two, the 
mass monopole and spin dipole. When these multipole 
moments are vifritten as a set of mass multipole moments 
Ml, which are nonzero for even values of /, and a set 
of current multipole moments Si, which are nonzero for 
odd values of I, the no-hair theorem is captured by the 
relation H 



Ml + iSi = M{iay 



(1) 



where a = J/M is the spin parameter. 

Thanks to the no-hair theorem, this property of black 
holes in general relativity naturally leads to the expec- 
tation that all astrophysical black holes are described by 
the Kerr metric. Astrophysical black holes, however, will 
be neither perfectly stationary, nor exist in perfect vac- 
uum. Other stars, electromagnetic fields, and other forms 
of matter like dust and dark matter, will induce pertur- 
bations away from the Kerr metric. Such perturbations 
will induce nonzero deviations from Eq. ([1} that could 
lead to a violation of the no-hair theorem. However, if 



one makes the implicit assumption that such perturba- 
tions will be so small to be essentially unobservable, then 
one can argue that astrophysical black holes must be de- 
scribed by the Kerr metric. This is the assumption I 
make in this paper. 

There exists observational evidence for the presence 
of horizons in astrophysical black holes (see discussion 
in, e.g., Ref. Q), but a proof of the validity of the no- 
hair theorem is still lacking. This is why a concrete 
effort has been brewing for the past few years to de- 
velop model-independent tests using electromagnetic and 
gravitational- wave observations to determine the precise 
strong-field nature of black holes. Instead of focusing on 
particular gravity theories and introducing modifications 
of the Einstein-Hilbert action, these tests are designed 
in a phenomenological approach that encompasses large 
classes of modified theories of gravity and which is able to 
test many different theories simultaneously. In this case, 
the underlying fundamental theory is usually unknown, 
and insight is hoped to be gained primarily through ob- 
servations. 

Such strong-field tests can be classified into two groups 
(see Ref. 28] for reviews on this topic): 

• Gravitational-wave tests using the gravitational 
waves generated by stellar-mass compact objects 
in tight extreme-mass ratio inspiral (EMRI) orbits 
around a supermassive black hole |8Hl9|: 

• Electromagnetic tests using the radiation emit- 
ted by accelerating particles in an accretion disk 
around a black hole |20l - l27| . 

Other weak-field tests of the no-hair theorem exist, such 



as those obtained from observing close stellar orbits 
around Sgr A* [23| and pulsar/black- hole binaries ,3Qj , 
but these do not probe the near-horizon, strong-field na- 
ture of black holes. While extensive searches of pulsars 
orbiting around black holes are ongoing (e.g., Ref. p3l|). 
no pulsar/black-hole binary has been found so far. 

In contrast to weak-field tests of gravity, where a pa- 
rameterized post- Newtonian approach (e.g., [323) i^ suffi- 
cient, strong-field tests require a careful modeling of the 
spacetime by introducing a parametric deviation from 
the Kerr metric. Several such parametric frameworks 
have been constructed, within which possible observa- 
tional signatures of deviations from the Kerr metric can 
be explored (e.g., [HEHIsi^). The objects they 
describe have spacetimes that can deviate slightly to 
severely from the Kerr metric, and observables can be 
studied in terms of one or more free parameters. All of 
these metrics contain the Kerr metric as the special case, 
when the deviations are dialed to zero. 

The many proposed metrics in the literature can be 
divided into two sub-classes: those that are Ricci flat. 



Rn 



0, and those that are not. In the former case, 



the metric in the far-field satisfies the Laplace equation, 
and thus, when in asymptotically Cartesian and mass- 
centered coordinates, it can be expressed as a sum of mass 
and current multipole moments (see, e.g., R ef. 13611) . One 
can relate these moments to each other via nU, [m [l3 



Me + iSi == MiiaY + 5Mt + '^SSi , 



(2) 



where 5Mi and 5Se are mass and current multipole de- 
formations. For this class of metrics, the measurements 
of three or more multipole moments could then be used 
to test for deviations in the Kerr metric Q. 

When the metric is not Ricci flat, the above param- 
eterization of the metric in the far-field (as a sum over 
mass and current multipole moments that depend only 
on the (. harmonic number) is not valid. Such metrics 
generically arise from explicit or implicit modifications to 
the Einstein-Hilbert action, such as in dynamical Chern- 
Simons gravity I37H39J and in Einstein-Dilaton-Gauss- 
Bonnet gravity [40]. In these cases, it is not clear what 
the general structure of a modification of Eq. ^ would 
look like. 

If a gravitational-wave or electromagnetic measure- 
ment requires a Kerr deviation, then there are two pos- 
sible implications. One possibility is that the object ob- 
served is not an ideal astrophysical black hole. This could 
mean that either the black hole is not perfectly stationary 
or it is not in pure vacuum, indicating a failure of my as- 
sumption that these prosaic deviations are unobservable. 
More interestingly, the object may not be a black hole at 
all, but is instead a more exotic object [Il|, |41i| that per- 
haps violates cosmic censorship [l| . Another possibility is 
that four-dimensional general relativity is only approxi- 
mately valid in the strong-field regime, and thus, station- 
ary and vacuum black holes solutions are not described 
by the Kerr metric (see, e.g., Refs. [13, IH, 113, El ) • In 
this interpretation, if the object is otherwise known to 



possess a horizon, both the no-hair theorem and strong- 
field general relativity are invalid. 

In this paper, I analyze the properties of several para- 
metric deviations from the Kerr spacetime. I consider 
the quasi-Kerr metric of Glampedakis & Babak (QK 
121), the bumpy Kerr metric of Vigeland & Hughes (BK 
ll[[l6[), the metric proposed by Manko & Novikov (MN 
_34|). the modified Kerr metric of Johannsen & Psaltis 
(MK; [33) and the modified-gravity, bumpy Kerr met- 
ric of Vigeland, Yunes, & Stein (MGBK; [ll]). I aim 
to identify the manner in which their properties deviate 
from the special properties of the Kerr metric as a con- 
sequence of the no-hair theorem. 

First I point out that the QK, BK, and MGBK metrics 
have been constructed as linear deviations from the Kerr 
metric [HI, [H [S [11 , while the MN and MK metrics 
are non-linear deviations from the Kerr metric |34l . l35j . 
The MN metric is Ricci-flat, the QK metric is Ricci- fiat 
up to terms containing the quadrupole moment, and the 
BK metric is a vacuum solution of the linearized Ein- 
stein equations if the spin vanishes. The MK and MGBK 
metrics are not Ricci-fiat. On the other hand, the QK, 
BK, MN, and MK metrics are stationary and axisymmet- 
ric and are generally of Petrov type I, while the MGBK 
metric also admits an approximate Carter constant and 
is of approximate Petrov type D. All of these metrics are 
asymptotically flat, which I show explicitly in the case of 
the QK and MN metrics (see Refs. [H, [3| for the BK, 
MK, and MGBK metrics). I show, however, that the MN 
metric requires a coordinate transformation and a rescal- 
ing of the mass in order to reduce to the Newtonian limit 
in the non-relativistic regime. 

I proceed to investigate the nature of the central object 
in all five metrics. Using techniques from the numerical 
relativity literature, I calculate the location of event hori- 
zons. The horizon in each case can be expressed as a level 
surface of a particular scalar function (see Ref. l59]). For 
the cases I consider here, this function is governed by a 
differential equation of the horizon radius as a function 
of the polar angle, which I solve using both analytical 
and numerical methods. I also derive an approximate 
analytic expression of this equation for small perturba- 
tions away from the Kerr metric. I show that the QK 
and BK metrics harbor naked singularities as is the case 
of the MN metric |3j| , while the MGBK metric describes 



a black hole HSf^. I also show that the MK metric har- 
bors a naked singularity, which is located at the Killing 
horizon and can have either spherical or disjoint topol- 
ogy. I calculate expressions for the deviation parameter 
as a function of the spin that delineate the boundaries 
between the regions of the parameter space with these 
different topologies. 

Finally, I identify the regions of space where viola- 
tions of Lorentz symmetry or closed timelike curves ex- 
ist, which I find outside of the central objects of the 
QK, BK, and MN metrics. These regions are unphysical 
and have to be excised by introducing a cutoff radius, 
which can, therefore, limit the ability of these metrics to 



serve as a framework for observational tests of the no-hair 
theorem. They impact both EMRI observations in the 
gravitational-wave spectrum, as well as electromagnetic 
observations of accretion flows, since both depend sensi- 
tively on the behavior of the metric near the innermost 
stable circular orbit (ISCO; see discussion in Ref. [35|). 
Consequently, the QK, BK, and MN metrics can only 
be used for tests in the electromagnetic spectrum for 
moderately-spinning black holes so that the ISCO lies 
outside of the pathological regions. 

I show that the MK and MGBK metrics are free of 
such pathologies exterior to the naked singularity and 
the event horizon, respectively, and argue that these 
metrics are particularly suited for electromagnetic and 
gravitational-wave tests, respectively. In the case of the 
MK metric, a cutoff radius has to be introduced just out- 
side of the naked singularity. Since the ISCO lies outside 
of the naked singularity for all values of the spin and the 
deviation parameter [35| , the cutoff radius can always be 
chosen so that the ISCO still lies in the domain exterior 
to the cutoff. The existence of a Carter-like constant in 
the MGBK metric allows one to separate the geodesic 
equations, which facilitates the study of EMRIs in such 
spacetimes. The MGBK is, thus, a useful tool for de- 
veloping gravitational-wave models that can be used for 
tests of the no-hair theorem [T^, [H, [H . 

This paper is organized as follows: In Section II, I com- 
pile the explicit forms of the five metrics that I study in 
this paper. In Section III, I discuss their symmetries and 
show that they are asymptotically flat. I analyze in detail 
the presence of event horizons in the five metrics in Sec- 
tion IV and identify pathological regions in Section V. I 
formulate my conclusions and discuss astrophysical appli- 
cations in Section VI. Throughout, I use geometric units, 
where G = c = 1. 



II. PARAMETRIC DEVIATIONS FROM THE 
KERR METRIC 

In this section, I summarize the explicit form of the five 
spacetimes [ij, [la, [la, [SJ, [3a| that I use in this paper to 
investigate parametric deviations from the Kerr metric. 

My starting point is the Kerr metric g^^, which in 
Boyer-Lindquist coordinates takes the form (e.g. (43|) 



K 



2Mr 
2Mar sin^ 



A 



2 I z 

r + a 



199 — ^ y 

2Ma'^rsm^e 



sin^ 6 . 



where 



A; 



2Mr- 



r^ + a^ cos^ I 



(3) 



(4) 



The only parameters in the Kerr metric are the first two 
multipole moments, i.e., the mass monopole Mq = M 
and the spin dipole Si = J — aM of the black hole. All 
multipole moments of higher order are related to mass 
and spin by Eq. ([ij. 

Some of the metric deformations that I study in this 
paper have been designed as linear deviations of the Kerr 
metric, which are of the form 



dfj.u ~ Sau + C "■/il' : 



(5) 



where C is a small parameter that reminds us that the 
metric deformation C^/^f is supposed to be small relative 
to the Kerr metric g^^,. These include the QK [12], BK 
[H [Hi, and MGBK [H] metrics. Unless I state oth- 
erwise, I perform my analysis of these metrics to linear 
order in the parameter C, i.e., I expand all quantities that 
derive from a metric of the form given by Eq. ([5]) to ©(C)- 
Other parametric deviations from the Kerr metric need 
not be linear and can be more general functions of the 
deviation parameters. These include the MN [3J| and the 
MK 35] metrics. I will analyze the properties of these 
metrics without expanding in the parameter C unless I 
state it explicitly. 

For the study of astrophysical phenomena in the met- 
rics with linear deviations from the Kerr metric, however, 
it is sometimes useful to compute their properties to all 
orders in the parameter ^, i.e., without expanding the 
results of such computations to 0{C,). While an expan- 
sion in the small parameter C, can always be performed 
in analytic calculations, it is a lot more difficult and, in 
some cases, even impossible to enforce in other settings 
such as the ones involving magnetohydrodynamic simu- 
lations, which numerically solve the (non-linear) geodesic 
equations. For this purpose, I also study the existence 
of event horizons and pathological regions in the QK, 
BK, and MGBK metrics to all orders in the parameter ^. 
Note, however, that in this paper I only consider small 
values of the parameter C,. My results for the QK, BK, 
and MGBK metrics without expanding in the parameter 
C, may be altered if larger values of the parameter C, are 
considered. Similarly, it is instructive to study the MN 
and MK metrics also in the limit of small deviations, 
expanding these metrics to first order in the deviation 
parameter and treating the resulting metrics as pertur- 
bative. 



A. The Quasi-Kerr Metric 

The QK metric [12j derives from the Hartle-Thorne 
metric [44| . which was originally designed for slowly ro- 
tating neutron stars. The QK metric deviates from the 
Hartle-Thorne metric in that its quadrupole moment is 
corrected, i.e., it is not assumed to depend on mass and 
spin through Eq. ([1]). Glampedakis & Babak 12| and Jo- 
hannsen & Psaltis [201 calculated orbital frequencies in 
this spacetime, and Johannsen & Psaltis [20|,l22] analyzed 
the properties of this spacetime including the locations of 



the ISCO and the circular photon orbit, the gravitational 
lensing experienced by photons, as well as the dynami- 
cal frequencies of thin accretion disks around the central 
object. 

The QK metric [12| modifies the quadrupolc moment 
of the Kerr metric by the amount 



SQqk 



-eQKM\ 



(6) 



where the parameter eqK measures deviations from the 
Kerr metric. The full quadrupole moment is then 



Qqk = -M {a^ + e^kM') 



(7) 



In Boyer-Lindquist-like coordinates, i.e., in spherical- 
like coordinates that reduce to Boyer-Lindquist coordi- 
nates in the Kerr limit, the QK metric g^^ is given by 
Eq. (O with Cqk = Cqk and [12] 

h'^^ = (1 - 2M/r)-i [(1 - 3cos2 0) J^,{r)] , 
/i^; = (1 - 2AI/r) [(1 - 3cos2 9) ^i(r)] , 



QK 



[(l-3cos2 0)j-2(r)], 



1 



QK 9 • 2 , 

r^ sm I 



[(l-3cos2 0)j-2(r)], 



(8) 



and /iqk = 0. The functions Ti^2ir) are given in Ap- 
pendix A of Ref. [13 • The QK metric is a solution 
of the vacuum Einstein equations for spins that satisfy 
a/M < 1, provided Eqk 7^ 0. 



B. The Bumpy Kerr Metric 

The BK metric [ll|, ll6| modifies the Kerr space- 
time through small perturbations due to external 
stresses. Collins & Hughes jllj] defined mass perturba- 
tions by starting from the most general stationary and 
spherically-symmetric metric in Weyl form 

ds^ = ^e^^dt^ + e^^-^^{dp^ + dz^) + e-^^p'd<p\ (9) 

They defined ip = iPqA- -01 and 7 = 70 -I- 71, with (V'Oi 7o) 
equal to the Schwarzschild values and ('01, 7i) parametric 
deformations. Both V'l and 71 are proportional to the 
small parameter Cbk- They then required this metric 
to satisfy the Einstein equations to linear order in Cbk, 
thus obtaining differential equations for the deformation 
functions. 

This was generalized to a spinning spacetime |16j 
by applying a Newman- Janis rotation [45|. In Boyer- 
Lindquist-like coordinates, the BK metric is given by 



Eq. © with [3 



/.ff --2(1- 



2Mr\ , 
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2aMr sin^ 9 



/jBK 
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^i)| 
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AS 
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2(71- 
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(71 
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^i''AS(S-2Mr) 


-2^Al- 
\ 


2Mr\ 



X A sin^ 9. 



(10) 



The perturbations ■0i and 71 satisfy linearized Einstein 
equations, which can be solved through a multipolar de- 
composition. At lowest order {£ = 2), these functions are 
given by [16| 



rr'i' 



71" 



B2M^ 



1 



B2\ - 



TT d{r, 9, a)^ 
L{r,0,a) 



3L{r, 9, a)2 cos^ 9 



d{r, 9, a) 



-1 



C2o{i", a) + C22{r, a) cos^ 9 + C2i{r, a) cos'' 
d{r,9,af 



-- 1 

(11) 



where 



2Mr+{M^ + a?)cos? 



d(r,0,a) = V^ 
i(r, 9, a) = V(r - Mf + a? cos^ 6*, 
C2o(r,a) == 2(r - Mf - bAP{r - Mf + 3Af^ 
C22(r, a) ^ 5M^{r - Mf - 3M'^ + a^[4(r - Mf - 5M^ 



C24(r,a) =a^{2a^ + 5M'^). 



(12) 



The strength of the perturbation to the Kerr metric at 
this order is determined by the parameter ^bk = ^2- 
Note the absolute value signs in the function d{r,9,a), 
which are missing from the corresponding expressions in 
Ref. [ia], Eq. (5.6). This function is the translation into 
Boyer-Lindquist coordinates of the Weyl-sector function 
cosh u cos^ V + sinh u sin v. It is positive definite in the 
Weyl sector, and should be positive definite in Boyer- 
Lindquist coordinates as well; '0i can become imaginary 
otherwise. 

Vigeland & Hughes [la| analyzed orbits and orbital 
frequencies in this spacetime, and Vigeland [13 showed 
that, at lowest order, the perturbation changes the mass 
quadrupole, which is given by [17| 



Q^^^-Mia^ + ^B2M^- 



so that 



SQbk 



1 



-B2M' 



(13) 



(14) 



This equation suggests a formal relationship between 
the deformation parameter Eqk of the QK metric and the 
parameter B2 of the BK metric via the relation 



:B, 



0.63B2. 



(15) 



This mapping is somewhat misleading, however, as one 
might be tempted to conclude that the QK and BK pa- 
rameterizations are equivalent when in fact they are not. 
The respective perturbations of the Kerr metric in the 
QK and BK metrics are different due to the different 
functional forms of the corrections h^^, and 6^^ given by 
Eqs. ([5]) and ([TU)) . respectively. In addition, for nonzero 
deviations from the Kerr metric, the QK metric is a 
solution of the vacuum Einstein equations only up to 
quadratic order in spin, while the BK metric is a solution 
of the vacuum Einstein equations, linearized in Cbk, to all 
multipole orders /. 



Following Ref . []J| , I define 



Qt.- 



M2 - M^ 



(17) 



which is a dimensionless parameter that measures the 
MN deviation from the Kerr quadrupole moment M^. 
The first few non- vanishing multipole moments are given 
by IM] (see, however, Ref. [Ij] and references therein) 



Mo = M, Si = aM, 

M2 = -M(a2 + gMNM2), 



(18) 



Expanding the MN metric to first order in the param- 
eter qmn, I obtain the linearized MN metric and identify 
Cmn = 9MN- Due to the lengthy form of this metric, I do 
not write it here explicitly. 



The Ricci-flat metric proposed by Manko &c 
Novikov 



The MN metric IM] is a non-linear superposition of 
the Kerr metric with a static vacuum field that general- 
izes the former to a Ricci-flat spacetime with arbitrary 
mass and current multipole moments. In its original 



form 3j|, this metric harbors a naked singularity 34 1. 
Gair et al. [1J| and Brink [Ij] analyzed properties of or- 
bits in the MN spacetime and found regions near the 
central singularity where geodesic motion becomes er- 
godic. They identified domains containing closed time- 
like curves around the origin which violate causality. In 
addition, Gair et al. [Till computed the location of the 
ISCO and the orbital frequencies. Berti et al. [4a] an- 
alyzed the MN metric in the context of rotating neu- 
tron stars, and Bambi & Barausse [2^, |43] and Bambi 
& Lukes-Gerakopoulos [4g] studied accretion disks and 
their thermal emission and potential gravitational-wave 
signatures in an MN background, respectively. Ergodic 
orbits in general non-Kerr spacetimes were analyzed by 
Refs. 0,l4|. 

In this paper, I will use a subclass of the MN metrics 
that describes electrically neutral compact objects and 
that is given by the line element of [34] transformed to 
Boyer-Lindquist-like coordinates 



dstt, = -/mn dt'^ + 2/mn w dt d(j) 



o2r 



(r - M)2 - (Af 2 _ q2) pog2 



Jmn 

^2r 



-dr^ 



+ ^[{r~ Mf - (M2 - a") cos2 9] dO^ 

JMN 

+ f-^ {AsinH - fl^u:') d4>^ 



(16) 



with /mn = fmNir,9), Lo = uj{r,6). The quantities T, /^, 
and uj are given in Appendix A. 



D. The modified Kerr metric proposed by 
Johannsen & Psaltis 



The metric proposed by Johannsen & Psaltis [35| 
contains a set of free parameters which introduce non- 
linear deviations from the Kerr metric. Johannsen & 
Psaltis [391 introduced polynomial corrections to the 
{t, t) and (r, r) elements of the Schwarzschild metric and 
transformed this ansatz into a Kerr-like metric via the 
Newman- Janis algorithm [4^. This metric is not Ricci- 
flat. 

The properties of iron lines, quasi-periodic variabil- 
ity, continuum spectra, and X-ray polarization from ac- 
cretion disks in the MK metric have been analyzed in 
Refs. [2J, [23, [5^. Moreover, the topology of this metric 
as well as its implications for the properties of black-hole 
jets have been studied in Refs. [S^, [51[ and [S^, respec- 
tively. Other properties of this metric were analyzed in 

Refs. Hi,!!!]. 

In Boyer-Lindquist-like coordinates, the metric is given 
by the expression 
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where 



h{r,9) = J2 (e2fe + e2/c 

k=l ^ 



Mr 



M2 



(19) 



(20) 



with free parameters e^. I will use this metric with only 
one non-zero parameter, so that the function h(r,9) re- 
duces to 



E. The Modifled Gravity Bumpy Kerr Metric 



h{r,9)^e3 






(21) 



Linearizing the MK metric to first order in the param- 
eter Cmk = £3,1 obtain the correction to the Kerr metric 
given by the expressions 
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The MGBK metric, proposed by Vigeland, Yunes, & 
Stein [18], deforms the Kerr metric through certain bump 
functions, such that the resulting metric possesses four 
constants of the motion. Such a metric is also not Ricci 
flat. Vigeland, Yunes, & Stein [ll] analyzed orbits in 
this spacetime and showed that speciflc choices of the 
bump functions reproduce all known modified-gravity 
black hole solutions known to date. Approximate EMRI 
waveforms in this metric were constructed in Ref. [19| . 

In this paper, I use the class of the MGBK metrics 
studied in Ref. [19|, where certain simplifications are 
made to guarantee certain properties (see Section IIB 
in Ref. llB]). With this at hand, the non- vanishing com- 
ponents of the MGBK metric in Boyer-Lindquist-like co- 
ordinates are as in Eq. ([S]) with 
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r 



where 



r^ - (? cos^ ( 



(r 



2 ' a2)2-a2Asin^6i. 



(24) 
(25) 



and Pj are polynomials in (r, cos ^) , given in the appendix 
of Ref. [1^ (I adopt here the deformed Kerr parameteri- 
zation of Ref. ^). 

The bumpy functions 7^ = 7i(r) depend on radius, and 
I parameterize this dependence via [l9| 



00 
7A = X! ^^^' 

n=0 



73 



- CXD 

- E T3." 



n=0 



(26) 



where A is 1 or 4 and (7i,n,73,Ti:74,n) are constants that 
control the magnitude of the deformations. Additional 



73,0 



simplifications [19] allow us to set 71^0 = 7i,i 

74,0 = 74,1 = 73,2 = 0. 

The metric components h^^^^"^ and h^?^^ are fully de- 
termined once I solve the differential equation in Eq. ((23|) 
for htcf,. Doing so in a far-field expansion, I find 



N 
n^2 



■t(f),n 



io) 



E"/2 



(27) 



where the coefficients /it0,n are given in Ref. [19]. Notice 
that the expressions presented here defer slightly from 
that of Ref. [19[ , as I use here a dimensional Kerr spin pa- 
rameter a and define S and p differently than in Ref. (19| . 
I likewise write the bump functions as 7^ instead of 7^ 
and the polynomials as Pj instead of Pj to avoid confu- 
sion with the bump functions in the BK metric and the 



Legendre polynomials that occur in the MN metric (see 
Appendix A). 

In this paper, I study the lowest-order perturbations 
and only allow the coefficients 71^2, 73,1, 73,3, and 74,2 
to be nonzero. This choice corresponds to setting TV = 2 
in Eq. ^. 



III. SYMMETRIES AND ASYMPTOTIC 
FLATNESS 



which the r -^ 00 limit is taken (see Ref. [58J for fur- 
ther details). However, this definition will suffice for my 
purposes in this paper. 

For all of the metrics described in Section [TTl the Kerr 
part is clearly asymptotically flat; I are thus left with the 
task of showing that the deviations of these metrics from 
the Kerr metric do not spoil the asymptotic flatness of 
the background. The asymptotic flatness of the BK, MK, 
and MGBK metrics has already been shown [I3, [Sa] ■ I 
now turn to the QK and MN metrics. 



By construction, all the metrics described in Section HIl 
admit two Killing vectors corresponding to stationarity 
and axisymmetry [III, IS Mj ISJ: IsJ • This implies that 
these spacetimes possess an associated conserved energy 
and conserved (z-component of the) angular momentum. 
The MGBK metric possesses in addition an approxi- 
mately conserved fourth quantity (a Carter-like constant) 
associated with the existence of an approximate Killing 
tensor (approximate in the sense that it satisfies Killing's 
equation to linear order in the deformation parameters) . 
Therefore, the MGBK metric is of approximate Petrov 
type D, while the other four metrics (QK, MN, BK, and 
MK) are, in general, of Petrov type I. 

Some metrics considered here satisfy the vacuum Ein- 
stein equations, while others do not. The MN metric is 
Ricci-flat, and thus, it is a solution of the vacuum Ein- 
stein equations. The QK metric is Ricci-flat only up to 
second order in spin and flrst order in the perturbation 
parameter, i.e., neglecting terms of ©(eqKo), 0{a?') and 
©(eqk). The BK metric is a vacuum solution only of 
the linearized Einstein equations in the limit a — )■ 0, i.e., 
the Ricci tensor contains terms of 0{aB2) [ia|. The MK 
and MGBK metrics do not satisfy the vacuum Einstein 
equations and are not Ricci-flat. 

In order to make meaningful predictions of observ- 
ables, spacetimes of black holes in isolation must be 
asymptotically flat, i.e., there must exist a coordinate 
system (a;°,a;^,a;^,a;^) such that all metric components 
in these coordinates behave as g^i, ~ rj^^, + 0{l/r) 
as r — ^ 00 in either spatial or null directions, where 
rj^i, — diag(— 1, 1, 1, 1) is the Minkowski metric and where 
r is the Euclidean norm of the spatial coordinates |55j . 
In terms of the line element in Boyer-Lindquist-like coor- 
dinates, the subleading terms must fall-off as (e.g., |56j'l: 



ds' 



- 


r ^ ' 


dt^ 




4: Ala / TX 

sin^ e + O (r-^) 


dtd(l) + 


dr^ + r^ {de^ + sin^ 9d(j)^) 





l + Oi 



(28) 



Asymptotically flat spacetimes with a slower fall-off can- 
not be stationary in general relativity [57|. The above 
definition of asymptotic fiatness is not precise, due to is- 
sues with coordinate invariance and the precise way in 



1. QK Metric 

As an example of the QK metric, I consider the (i, t) 
component, which has the form 



hft'^ ex -^{2M^ + AM^r - 9Mr^ + Sr^) 



- 3(r - 2M)2 In 



r - 2M 



(29) 



One can expand the logarithm in r 3> TV/ to show that 
hf^ ex — (16/5)M^/r^, which clearly remains asymptot- 
ically fiat. A similar argument holds for the other com- 
ponents of the QK metric. 



2. MN Metric 

The MN metric should be considered separately, as 
here its asymptotic fiatness is not as obvious. To see its 
structure more clearly, I perform the coordinate trans- 
formation 



r' = exp 



4g^ 



(1_ ^2)3/2 



exp 



4q^ 



(1 - x2)3/2 



(30) 

with (i, 9) unchanged and x = ci/M. This transformation 
ensures that the MN metric reduces to the Minkowski 
spacetime at radial infinity. Transforming the metric, 
expanding its elements in M/r ^ 1 and linearizing them 
in gMN <C 1, the metric perturbations become (dropping 
primes) 



7 MN 



LMN 
"•40 



T^MN 

'Vr 



,M 1 



.2^3/2 



(1-x 

32xgMNM2 sin^ 9 1 



-+0{r-'), 
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SOmnM 1 
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X^ sin^ 
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(1 - X^f 

mnx'^' sin' 9 
(I-X'f^ 



o{ 



,-3^ 



+ 0(r-i). 



Oir- 



(31) 
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The parameter M has to be rescaled in order for the 
MN metric to describe the correct Newtonian hmit in 
the non-relativistic regime: 



M' = M exp 



4g^ 



(1_ ^2)3/2 



(32) 



This then ehminates all components of the metric per- 
turbation to relative 0{r~^), with the deformed met- 
ric becoming equal to the Kerr one with mass M' and 
spin a' = xM', and shows that the MN metric is indeed 
asymptotically flat. An investigation of a rescaling of 
the fuU metric and its parameters at O {M^ / r^ , ql^^^) is 
beyond the scope of my analysis. From here on, I will 
use the MN metric in the form given by Eq. (IT5|) after 
applying the coordinate transformation in Eqs. pop and 
the rescaling in Eq. ([5^ . 



IV. THE EVENT HORIZON 

The event horizon of a black hole delineates the re- 
gion of spacetime which cannot communicate with dis- 
tant observers (its interior) from the region which can 
communicate. In this section, I describe the calculation 
of the location of event horizons in stationary, axisym- 
metric, asymptotically flat metrics such as the ones listed 
in Section II, using techniques developed in numerical rel- 
ativity. I proceed to compute the location of the event 
horizons of these metrics should they exist. 

The event horizon is a null surface, generated by null 
geodesies ("generators") that are trapped within that 
surface. The normal to a null surface, n'^, is itself null: 
n^Ti^ = 0. I can take this surface to be the level surface 
of a scalar function f{x°'), in which case the normal is 
simply n'^ = V^/ — d'^f. The event horizon is then de- 
fined by the condition (see, e.g., Ref. [531 ^^^ ^ detailed 
discussion) 



J^'id^Did.f) = 0. 



(33) 



Choosing a time coordinate, this becomes a quadratic 
equation for dtf which can be solved to show how the 
horizon evolves given an initial condition. This equation 
can also be evolved backwards given some final condi- 
tion. Very powerful codes have been developed in the 
past decade which solve the "master" horizon equation, 
Eq. (|33l) . in dynamical spacetimes. As long as the solu- 
tion settles down to the Kerr metric at late times, these 
codes can find the level surface / that describes the hori- 
zon quite robustly ;f59J. 

One can distinguish between coordinate singularities, 
coordinate locations where the line element diverges, and 
true singularities by evaluating curvature scalars. I will 
here distinguish between these two cases by computing 
the Kretschmann scalar. 



K = R^p^sR"^''', 



(34) 



The metrics described in Section II are all parameter- 
ized by spherical-like coordinates (t, r, 9, </>) and are sta- 
tionary and axisymmetric. In this case, the function / 
which characterizes the horizon can then only depend on 
the coordinates r and 9, and I have 

g'^idrf)' + 2g^\drf){def) + g'^deff - 0. (35) 



Eq. ([55]) defines event horizons in the spacetimes that 
I study. Note, however, that the existence of a solution 
of this equation is only necessary for the presence of an 
event horizon, but generally not sufficient, because solu- 
tions need not form a closed surface or can be singular, 
i.e., the Kretschmann scalar can diverge at this location. 
In these cases, a solution of Eq. (l35t is simply a null 
surface, but not an event horizon. In this paper, I will 
distinguish between the two as appropriate. 

In the following, I will choose special forms for the level 
surface function /, which further simplify my analysis. I 
motivate this special form by first examining the Kerr 
spacetime, and then generalizing this to the spacetimes I 
consider here. 



A. Kerr Black Holes 

A useful first guess for the scalar function / is the radial 
coordinate r. Since / is only defined up to a constant 
[only its derivatives enter Eq. (l33l) ]. it is useful to subtract 
a constant so that / = on the horizon. Let us then set 



/ 



H, 



(36) 



where H is the (initially unknown) location of the space- 
time's event horizon. Level surfaces of this function / 
are a sequence of nested coordinate spheres, and Eq. ((55)) 
simplifies to 



5-(iJ) = 0, 



(37) 



which follows from choosing f — 0, or simply r ^ H , to 
define the horizon. 

It is sometimes erroneously thought that this condition 
generically describes event horizons. In fact, this is only 
true if surfaces of constant r have uniform causal struc- 
ture, i.e., if constant r surfaces are everywhere spacelike, 
null, or timelike (assuming that they are closed). If this 
is not the case, which depends on the underlying coor- 
dinate system, then this condition will give the wrong 
solution. If constant r surfaces do have uniform causal 
structure, then I find that level surfaces of / are space- 
like at large radius, timelike at small radius, and null at 
r = H. The classic black hole solutions of general relativ- 
ity are of this type, at least in the standard coordinates 
used to describe them. The solution in this case is given 
by 



H = Hk ^ r+ = M + Vm^ - a2 



(38) 



where R^a^g is the Riemann tensor. 



(ignoring the possibility of charge). Note that the Kerr 
horizon radius is typically denoted by r+, while H is 



often used in the numerical relativity literature to denote 
the horizon radius, which can vary with both time and 
position. I will use both notations in this section. 

For any stationary and asymptotically flat spacetime, 
the event horizon is also a Killing horizon for some Killing 
vector x'' if Hawking's rigidity theorem [3| or perhaps an 
appropriate generalization thereof (see, e.g., Ref. [60J) 
can be applied. Then, this Killing vector can be written 
as x^ = t^ + f2 (j)'^ for some constant 0, where t^ and 
(j)^ are the temporal and azimuthal Killing vectors of the 
spacetime (see, e.g., Ref. [61|). Requiring that this vector 
be null forces the condition 



.900 I ^ 
which I can use to solve for the constant 



g<p<p 94,4, ) 



.900 V -900 500 



(39) 



(40) 



As one approaches the event horizon, the angular velocity 
r2± must approach a constant, r2± — )■ — 5t0/.g00. This is 
because this constant represents the angular velocity of 
zero-angular momentum observers at the horizon, which 
must be single- valued. The only way this can happen is 
if (see, e.g., Ref. [H) 



5t0 - 5ttff00 — ■ 



(41) 



The radius at which Eq. (|1T|) is satisfied defines the 
Killing horizon for the spacetime, since there x^ is null. 
Eq. (|1T|) is then equivalent to Eq. ([57|. In general, how- 
ever, the Killing and event horizons of a stationary space- 
time are distinct (see, also, Ref. ^62j). 



B. Existence Conditions for Horizons of Modified 
Black Holes 

In the spacetimes I consider, the metrics are more com- 
plicated than the Kerr metric, and the f = r — H ansatz 
is not sufficient. I instead let the horizon radius be a 
function of 9: 



f = r^H{e). 



(42) 



This form is sufficiently general to describe any horizon 
for which there is a unique horizon radius for any given 
angle 9; see Ref. [531 fo'^ discussion. This ansatz may 
not be adequate for extreme deformations or for horizons 
with unusual topology (e.g., toroidal or disjoint horizons; 
see Fig. 5 of Ref. ,63] for an example). It will, however, 
suffice for my discussion in this paper as I will show in 
this section. 

With this in mind, the horizon is then defined by the 
condition r ~ H{9), where, using Eqs. (1551) and (H^ . 
H{9) is defined by the condition 



2r I '4 



dH 
~dB 



= 0. 



(43) 



The metric components which appear here are each 
functions of r and 9. I evaluate them at r = H{d), and 
so Eq. (|43|) is an ordinary differential equation for the 
horizon radius H{9). A formulation of this kind is used 
in most numerical horizon finders, modulo some small 
modifications to handle horizons of unusual topology (see 
Refs. [59i, 163!] for detailed discussion). 

I next examine how to solve this equation for the par- 
ticular cases that I study. In all of these cases, I work in 
coordinates such that g^^ = 0, so I drop the cross term 
in Eq. (|^5| in what follows. 



C. Linearly Deviating Spacetimes 

The QK, BK, and MGBK spacetimes are only specified 
as linear deviations from the Kerr metric, see Eq. ([5]). 
In order to analyze the existence of an event horizon in 
these metrics to linear order in the respective deviation 
parameters, I expand the function H{9) as 



H{9)=Hk+CSH{9), 



(44) 



where I introduced the Kerr horizon radius Hk here so I 
automatically find the Kerr solution for C = 0. Using 



9' 



,f.LU LiU 



Ch^'", 



(45) 



Eq. ((32]) becomes 



97 - Ch^'^ + i9'J - Ch'') (c^J 

Truncating at linear order and using h^^ 
further simplify this equation to 



(46) 



9k 9k "-r 



(47) 



which is equivalent to Eq. ([ST)) at 0{C,). I now examine 
what this equation implies for the QK, BK, and MGBK 
metrics. 



1. QK Metric 

In the QK metric, the element g^Tj^ is given by the 
expression 



..rr 



4+ Sid -3™=, 



2 (3r3 - %Mr^ + iM^r + 2M^ 



-3r 



2/2 
'r 



2M^] In 



r-2M 



(48) 



If £qk 7^ 0, Eq. (|47|) does not have a solution for all 



values of the angle 9. Solutions only exist as long as 
r > 2M; at radius r = 2M the logarithm in Eq. pS)) 
diverges. Evaluating the Kretschmann scalar K given 
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by Eq. (jMl) and expanding K to ©(cqk), I find that K 
diverges at the radius r — 2M, which I, thus, identify as 
a singularity. The solution of Eq. (j47|) . therefore, forms 
a null surface and not an event horizon, and the object 
is a naked singularity. 

For positive values of the parameter eqK, the null sur- 
face (where present) has a more oblate shape relative 
to the horizon of a Kerr black hole of equal spin, while 
the shape of the null surface is more prolate for negative 
values of the parameter Eqk. The location of the null 
surface in the equatorial plane increases with increasing 
values of the parameter Eqk as already found in Ref . [20] . 
There, however, the null surface was erroneously identi- 
fied as the Killing horizon determined by the condition 
(HTI) leading to a slight difference in the location of this 
surface. In Fig. [U I plot, for illustrative purposes, the 
QK null surface for \a\ = 0.3M and several values of the 
parameter Eqk in the xz— plane, where x = vr^ + o^sin^ 
and z = r cos 9. 



2. BK Metric 

If I look at the BK metric in Eq. (fTO|) and the defini- 
tions of the perturbation functions '0i and 71 in Eq. (jlll) , 
I see that the metric becomes singular at r+ given by 
Eq. dill), as weh as where E- 2Mr = and d(r, 6*, a) = 0. 
The second condition occurs when r is given by 



ri,± =M± VM2 



(49) 



where I note that ri^-|- > r"!,-, and ri^^ coincides with 
the location of the Kerr ergosphere. The third condition 
occurs when r is given by 



r2.,± 



M ± v^Af 2 sin^ I 



(50) 



where clearly r2^+ > 7^2,- • The Kretschmann scalar K 
given by Eq. ([M)) expanded to C'(i?2) diverges at the 
singularities ri_± and r2.±, which are, therefore, real sin- 
gularities. 

To leading order in the perturbation, the element g^^ 
is 



5bk 



^[l-2(7i-Vi)]- 



(51) 
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FIG. 1. Null surface of the central object in the QK metric 
with a spin of \a\ = 0.3M for several values of the parameter 
eqk. If eqK 7^ 0, the null surface is bound by the singularity 
located at radius r = 2Af , and the central object is a naked 
singularity. 



In their construction of the QK metric, Glampedakis 
& Babak ^2^ limited the validity of their metric to ex- 
clude the central region where r < 2Af due to the sin- 
gularity located at radius r — 2M . A corresponding 
cutoff radius ^cutoff (a, egx) > 2A'f as a function of the 
spin and deviation parameter was defined heuristically 
in Ref. |20| denoting an inner boundary of the region 
where the QK metric provides a consistent description of 
spacetime without pathologies. 



Consequently, as long as the bump functions 71 and V'l 
are regular at A = 0, then this last equality signals the 
location of a null surface. If, for some reason, the bump 
functions are singular here, then the null surface location 
would be modified from its Kerr value. Due to the pres- 
ence of the singularities, the null surface cannot be an 
event horizon. In Fig. [5] I plot the location of the null 
surface and the singularities in the xz-plane for a — 0, 
\a\ = 0.5Af , and \a\ — M. Since a appears in the expres- 
sions for ri ± and T2^± as a? , the locations of the singu- 
larities are the same for positive and negative values of 
the spin. For non-spinning black holes, the singularity 
ri^-i- coincides with the BK null surface. The singularity 
r-2,+ is hidden inside this surface at angles < 6* < 7r/2 
and coincides with this surface in the equatorial plane 
[9 — 7r/2). As \a\ increases, the singularity ri_+ moves 
outside of this surface. For small but nonzero values of 
the spin |a|, the null surface is closed if -82 > and ter- 
minates at the singularity r2.+ near the equatorial plane 
if -82 < (see Figure [8] for an example). This equatorial 
"hole" increases as |a| increases. 

The existence of naked singularities is consistent with 
the picture of the burupy black holes originally described 
by Collins & Hughes [il[ , in which naked singularities of 
the Curzon type were explicitly introduced to chan ge a 
spacetime's multipolar structure. Collins & Hughes [11| 
called attention to this strong-field naked singularity, 
cautioning that one could find odd results by using their 
spacetime to study very strong field structures. The 
same cautionary note clearly applies to the Vigeland & 
Hughes [la] bumpy black hole as well. Though useful for 
studying phenomena sensitive to orbits with separation 
r > (a few) x ri._|_, the naked singularities that are intro- 
duced will surely have an adverse impact for orbits with 



11 



S 






0,5 1 L5 2 

X (M) 



FIG. 2. Location of the BK null surface and the singularities in the a;z-plane for a = 0, \a\ = 0.5M, and \a\ = M. For a = 0, 
the singularity ri,+ coincides with the null surface, and the singularities r2,+ and r2,- are located inside the BK null surface. 
At 9 = 7r/2, the null surface intersects with the singularity r2,+ . As \a\ increases, the singularity ri,+ moves outside of the BK 
null surface. For small but positive values of the spin |a|, the null surface is closed if B2 > 0. Otherwise, it terminates at the 
singularity r2,+ . For all values of the deviation parameter B2 ^ the object constitutes a naked singularity. 



3. MGBK Metric 



For the MGBK metric, to leading order in the pertur- 
bation, the element 5mqbk i^ 



5MGBK = ^(l+7i), 



(52) 



where 71 is given by Eq. (PB)) . Similar to the case of 
the BK metric, the null surface is, then, determined by 
either the condition 71 = —1 or by the condition A = 0, 
where the latter condition holds if and only if 71 remains 
regular at that radial location. In particular, if 71 „ > 
for n > 2, the null surface coincides with the Kerr event 
horizon r^^J^^ — r_|.. Moreover, all components of the 
metric and the inverse metric become singular only at 
r = r+ or at r = 0. However, by direct evaluation of 
the Kretschmann scalar, it is clear that r = r_|_ is the 
location of a coordinate singularity, and not an essential 
one (the latter remains at the Kerr ring singularity) . The 
null surface is, therefore, an event horizon. These results 
were also found in Section II C of Ref. fig]. In other limits 
of the MGBK metric (e.g., Refs. ^I6|l), the location 
of the horizon is different from the location of the Kerr 
event horizon, as I demonstrate explicitly in Appendix B. 



D. Spacetimes with Non-linear Deviations 

The MN and MK spacetimes are non-linear functions 
of the respective deviation parameters, and the full hori- 
zon equation (jlH]) has to be solved in order to study the 
existence of event horizons in these spacetimes. Likewise, 
one may be interested in the full solution of this equa- 
tion to all orders in the parameter C, for the QK, BK, and 
MGBK metrics. In this case, one no longer determines 
the location of a null surface or an event horizon using 



Eq. (|47p , which is the version of Eq. (P5)) that has been 
linearized in the deviation parameter C,. In this section, 
I solve Eq. (|43p numerically for the horizon radius H{9) 
for all five spacetimes. In Appendix C I describe two 
such methods, a spectral and a finite difference method, 
which I used for my analysis. First, however, I explore 
Eq. p3)) analytically taking advantage of the symmetries 
of these spacetimes. 

Due to axisymmetry and reflection symmetry, the nor- 
mal to the horizon must be purely radial at the poles 
[9 = 0, tt) and in the equatorial plane {9 = 7r/2). In 
these cases, the horizon equation, Eq. (H5)) . simplifies to 
Eq. ([37)) . which can be solved directly. 



1. QK Metric 

For the QK metric, I find that the null surface coin- 
cides with the singularity at r = 2M at the poles and 
in the equatorial plane. Since the metric contains terms 
ex ln[r/(r — 2M)], which diverge at these locations, the 
numerical methods described in Appendix C do not con- 
verge. However, I find another singularity located at 
the Killing horizon specified by Eq. (|4ip as long as it 
lies outside of the singularity at r = 2M. For values of 
the parameter Cqk > 0, the Killing horizon lies around 
the equatorial plane, while for values of the parameter 
Eqk < 0, the Killing horizon lies around the poles. 



2. BK Metric 

For the BK metric, I find that the null surface coincides 
with the singularity ri__|_ at the poles and that it does 
not pass through the equatorial plane if i?2 7^ 0. Due to 
the singularities, my numerical algorithms likewise did 
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not converge, but I suspect that the nuU surface, where 
present, hes inside of the singularity ri _(. as it does if the 
BK metric is treated perturbatively (at least for |a| > 0). 



3. MGBK Metric 

For the MGBK metric, the event horizon, again, coin- 
cides with the Kerr event horizon at the poles and in the 
equatorial plane. Numerically, I find that it is identical 
to the Kerr event horizon at all angles and that the 
Kretschmann scalar remains finite there. 



4. MN Metric 

In prolate spheroidal coordinates (t,x,y, (/>), the 
MN metric harbors a naked singularity [34|. In 
Boyer-Lindquist-like coordinates I find numerically that 
Eq. (|37| has no solution in the equatorial plane for any 
value of the deviation parameter ^mn ^ 0. I confirmed 
the absence of a horizon with the finite difference method, 
which did not converge if ^mn j^ 0. 

However, by evaluating the metric elements, I find that 
the MN metric becomes singular at the location of the 
Kerr event horizon r+, at which the elements g^^ and 
g^^ either diverge or vanish depending on the sign of the 
parameter qmn- Due to the lengthy form of the MN met- 
ric, I did not calculate its Kretschmann scalar to identify 
the nature of this singularity. For radii r < r-^-, the MN 
metric can become imaginary (see Figure [lU] for an ex- 
ample), and I, therefore, suspect the presence of a true 
singularity. My results remain if I linearize the MN met- 
ric to 0((7mn)- 



5. MK Metric 

The relevant components of the MK metric are given 
by the expressions 



9mk 



9mk 



A + a^sm'^0h{r,e) 



E[l + h{r,e)] 



1 



(53) 
(54) 



where the function h{r,9) is given by Eq. ((2T|) . First, I 
examine the metric element gJ^K ^^ the poles and in the 
equatorial plane, again taking advantage of axisymmetry 
and reflection symmetry, li 9 ^ or 9 ~ n, g^K ^ 
A. This equation has a root at r = Hk — ?"+ unless 
£3 — —Ar^/M, in which case the denominator likewise 
vanishes and gJ^K = 1/2. Therefore, the event horizon or 
null surface (if present) must pass through the z-axis at 
the Kerr horizon radius. 

I explore next how any null surface must behave as I 
move in 9 away from the z-axis. I examine H{9) for small 



angles, 69 <^1, putting 



H{69) = Hk + > 



dH 
~dB 



q2 j2 



d'H 



e=Q 



2 rf6i2 



(55) 



61=0 



Note that a similar expansion describes the behavior of H 
near 9 — tt; thanks to reflection symmetry, it is sufficient 
to focus on 6 = 0. 

I insert r = H{59) into Eq. (|43l) . evaluate all metric 
functions at 6* = 59, and expand in 59. To 0{59), I find 







dH 



?=o 



{dH/d9)\g^o 59 / VAf 2 - a2 



2MHk 



M \ H 



K 



Me3/4 



{d''H/d9^)\e=o [{dH/d9)\e=aY 



Hk 2MHk 

From this expression, I conclude that 



dH 
~dB 



= 



(56) 



(57) 



?=o 



in agreement with axisymmetry. 

At second order in 59, I find a quadratic equation for 
d'^H/d9'^: 



=a^Me3 + AiMHK - a^) 



d^H 



-{AHk + Me-i) 



d^H 



d9^ 



d9^ 



*=o 



This equation only has real solutions as long as 
a^MtziAHx + Mes) < "^{MHk - a^f . 



(58) 



(59) 



Otherwise, d^H/d9^ is imaginary, and a null surface does 
not exist. Eq. (|59|) implies an upper and a lower bound 
on the deviation parameter 63, given by the expressions 



niin— polo 



< £3 <e 



max— polo 



where 



max-polc _ ^V ^IVITIK - a^ "^Hk 



= 2 



niin — pole rj 



V2MHk 

|o| M ' 

V2MHk - a^ 2Hk 
\a\ M~' 



(60) 

(61) 
(62) 



In the equatorial plane, the condition given by Eq. ([37]) 
reduces to 



r^ - 2Mr 



,A-f3 



(63) 



There are no closed-form solutions to this quintic equa- 
tion, but it is simple enough to find a numerical solution. 
Doing so, I find generically that there is a maximum pos- 
itive value of £3 for which Eq. (j63p has a real solution. To 
understand this £3 threshold value, consider the shape of 



5mk( 



tt/2) as a function of r. For modest positive 
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values of €3, it has two roots which tend to be close to 
r± = M i: \JM'^ — a?- . There is also a minimum at 



TT. > (64) 




which is roughly halfway between these roots. As £3 is 
increased, the roots move toward the minimum, with all 
three points converging at a value g™'^'^""'^ given by 



TT I 

Solving this equation yields 
1 



/o max— cq\ p, 

= 7r/2,e3 ') = 0. 



(65) 



max— cq 



3125(a/M)2 



1024 U + yi6 - 15(a/M)- 



160(a/M)2 (40 + 7v/l6-15(a/A/)- 



150(a/Af )■* ( 15 + VI6 - 15(a/M)2 



(66) 



Therefore, a null surface can only exist for values of the 
deviation parameter £3 < g™'*'^~°i given by Eq. (1661) . The 
bound g™^''^°i^ however, is weaker than the condition 



max— cq 



> £ 



max— pole 



given in Eq. (pO]) . since e 

Imagine now that a null surface does in fact pass 
through the equator, and consider how it behaves as I 
move bd away from the equatorial plane. Let i?oq be the 
null surface radius at 9 = 7r/2 and put 



H{^ 12+56) = H,^+SO— 



<Ph 



e=iT/2 



2 de^ 



(67) 



?=7r/2 



I insert r = H{-k/2 + 89) into Eq. gSj) and put = 
7r/2 + 56. Note that in this case, I cannot do as much 
analytic exploration. Because -ffoq must itself be solved 
numerically, much of my analysis must likewise be nu- 
merical. 

I find that {dH/d6)\g^^/2 = as expected for 
all parameters examined and that a real solution for 
{d'^H/de^)\g^^/2 exists for ah values < £3 < e^'^'^-P°i'=. 
For £3 < and \a\ > 0.82M, however, there exists a part 
of the parameter space in the range g™™~P° ° < ^^ < 
within which there is no real solution for d^H/d6^. 
Therefore, a null surface cannot exist in this region ei- 
ther. 

In order to further analyze the nature of the central 
object in the MK metric, I calculate the location of the 
Killing horizon using Eq. (HTI) . Setting 6* = or 6* = tt in 
Eq. gl]), I find 



5r3r-(5.r)' 



ex A 



1+63 



"M^ 



(r^ 



(68) 



Therefore, the Killing horizon coincides with the Kerr 
event horizon (and hence with the null surface, if it exists) 
at the poles. For values of the parameter £3 < £3 ' '^° , 
where 



Kil-pol 



16 



(- 
3V3 VAf 



(69) 



Killing horizons in addition to the (outer and inner) 
Killing horizons emerge near the origin; one such hori- 



zon emerges if £3 = £ 



Kil-pol 



and two if £3 < e 



Kil-pol 



Eq. (|4T|) in the equatorial plane [6 — 7r/2) reduces to 



1-K£3^r I ( A-|-£3 



a^M-^ 



(70) 



where the second factor is identical to Eq. (|63|) , which de- 
termines the location of the null surface in the equatorial 
plane. This equation has two solutions if £3 < ^™^^~^^ 
[c.f., Eq. (|66|)]. which are identical to the locations of 
the null surfaces (i.e., the location of the inner and outer 
null surface if < £3 < f^'^^~°^ and of the null surface if 
£3 < 0, in which case there exists only one). If £3 — —8, 
the inner and outer Killing horizon coincide in the equa- 
torial plane, and if £3 < —8, these two horizons cross, 
i.e., the inner Killing horizon lies outside of the outer 
Killing horizon in and around the equatorial plane, while 
it lies inside near the poles. For values of the parameter 



£3 > £ 



max— cq 



the Killing horizon does not pass through 



the equatorial plane, and the Killing horizon is disjoint 
forming two spherical surfaces centered on the symmetry 
axis above and below the equatorial plane. 

Numerically, I find that a null surface exists if 0.1 > 
£3 > —8. In this region of the parameter space, the 
Killing horizon lies slightly outside of the null surface at 
the polar angles < 9 < 7r/2 and 7r/2 < 6* < tt if £3 ^ 0. 
Both surfaces coincide at all polar angles only if £3 = 0, 
in which case they form the Kerr event horizon as encap- 
sulated in Hawking's rigidity theorem [3| . Outside of this 
region, the existence of a null surface is uncertain due to 
increasing numerical error in the solution of Eq. (|43p . 

Evaluating the Kretschmann scalar of the MK metric, 
I find that it diverges at the Killing horizon at all angles 
< 6 < TT. Due to the polar coordinate singularity of 
the Boyer-Lindquist-like coordinates, it is unclear if the 
Kretschmann scalar likewise diverges at the poles of the 
Killing horizon. The Killing horizon exists for all values 
of the spin and the parameter £3 even if there is no null 
surface. Since the Killing horizon is singular if £3 7^ 0, 
the MK metric harbors a naked singularity located at the 
Killing horizon. 

In Figure |31 I plot the various regions of the parameter 
space of the MK metric. At a given value of the spin, the 
shapes of the Killing horizon (if it has spherical topology) 
and, if present, the null surface are more prolate than the 
Kerr event horizon for values of the parameter £3 > 0, 
while they are more oblate for values of the parameter 
£3 < (see Ref. [35|)- I plot illustrative examples of 
the topology transition of the Killing horizon across the 
boundary e^'^^~°'^(a) given by Eq. ()66|) in Figured 



The bound ej^^^ "'^('j) coincides with the bound found 
in Ref. |35| , which mistakenly assumed the null surface to 
be a Killing horizon and used Eq. (|4T|) instead of Eq. pS]) . 
As shown in Ref. [35|, for positive values of the spin 
this boundary also separates the regions of the parameter 
space where no ISCO exists [i.e., the blue region on the 
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FIG. 3. Region of the parameter space of the MK metric 
where a null surface exists. In this and the shaded regions, for 
values of the parameter £3 ^ 0, the central object is a naked 
singularity located at the Killing horizon, which is of spherical 
topology. If a null surface exists, the Killing horizon coincides 
with the null surface at the poles and in the equatorial plane 
and lies outside of the null surface otherwise. In the blue 
shaded regions, the Killing horizon is of disjoint topology. 
Due to numerical uncertainties in the determination of the 
location of the null surface, I find that a null surface may not 
exist in this region if es < —8 or if £3 ^0.1. The black dashed 
line corresponds to a Kerr black hole. 



right-hand side in Figure[3l where es > e^^^~'^°^(a) ^ a > 0] 
from the region where an ISCO exists (i.e., everywhere 
else in Figure [3]). If an ISCO exists, it is always located 
outside of the Killing horizon and, therefore, outside of 
any region containing singularities [35]. If an ISCO does 
not exist, the only singularity in the equatorial plane is 
the origin; I analyze the corresponding part of the pa- 
rameter space in detail in a separate paper [65[ . Ref. [51j 
used the condition given by Eq. (|37l) to focate the event 
horizon, which likewise did not take the proper angular 
dependence of the horizon into account. 

In this paper, I only discuss the various metrics in the 
spin range |o| < M, because the Kerr background itself 
is pathological outside of this range (see discussion in the 
next section). For completeness, however, I investigate 
the nature of the central object in the MK metric for 
values of the spin |a| > M observing another change in 
the topology of the Killing horizon. If \a\ > M, Eq. ([70)1 
has two solutions if £3 < and no solutions otherwise, 
while Eq. 



tion if ea = e 



has two solutions if £3 < eg ' ^° , one solu- 

Kil-pol 



and no solution otherwise. Solving 



Eq. ((4T|) at the remaining polar angles, I find that the 
inner and outer Killing horizons have spherical topol- 



Kil- 

3 



-pol 



ogy It £3 S £3"' ""' as already found in Ref. [35|- If 

£3 = £3 ' ~^° , the inner and outer Killing horizons coin- 
cide at the poles. If £3 > £3 ' ~^° , the Killing horizon 
has toroidal topology and is centered in the equatorial 
plane. In Figure [51 I plot several examples of the differ- 
ent topologies. Such changes in topology are similar to 
the ones reported in Ref. "5\]\. 

Finally, I investigate the nature of the central object 
in the MK metric when it is treated as a small pertur- 
bation of the Kerr metric in the sense of Eq. ([5]) . In this 
case, I can calculate the location of the null surface from 
Eq. (|47l) . which, in this case, turns out to be an event 
horizon. Note that the metric element h^^ in expression 
Eq. p2|) diverges as r — )• r+; however, h^^ = gUgHh^J^ 
is well behaved there, since g^ = A/E, canceling the 
factor A~^. Since this expression is still valid for all spin 
values \a\ < M, I expect the event horizon to be located 
at the radius 



rn = r+{l + A£3), 



(71) 



where r+ is the Kerr horizon given by Eq. ([38]) and where 
A is the amount, by which the horizon is modified rela- 
tive to the Kerr horizon. Inserting this expression into 
Eq. (|47|) and linearizing in the parameter £3, I obtain the 
equation 



A = -- 



e^a'^M^r^ sin 



2VM^ - a2 (2A/r+ - a^ sin^ 6) 
The event horizon is then located at the radius 



rn 



1 



£3a^APsin''6l 



2VM2 



■ (2A'f r_| 



(72) 



(73) 



The Kretschmann scalar of the linearized MK metric re- 
mains finite at all radii r > 0, and the central object is 
a black hole for both positive and negative values of the 
parameter £3. 



V. LORENTZ VIOLATIONS, CLOSED 

TIMELIKE CURVES, AND REGIONS OF 

VALIDITY 

Each of the metrics discussed in Section HIl may harbor 
regions of space where the Lorentzian symmetry is broken 
or which contain closed timelike curves. In Section IIII[ 
1 studied the asymptotic structure of the metrics at spa- 
tial infinity, but this does not guarantee that they retain 
their Lorentzian signatures close to the central objects. 
Since I already showed that the metrics have Lorentzian 
signature at spatial infinity since they are asymptotically 
flat, their signature must change sign close but outside 
their central objects. In that case, their determinants 
must vanish somewhere outside of their central objects, 
making them singular. If such regions exist, the metric 
can only describe physical processes outside of them. In 
the following, I analyze the properties of the five metrics 
in this regard. 
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FIG. 4. Different shapes of the Kilhng horizon in the MK metric for values of the spin \a\ — 0.9M and deviation parameter 
£3. Left: £3 = 0.31 < e™'"'~'"'(a), where e^'^^~'"^{a), a > 0, denotes the boundary between the regions of the parameter space 
with different topologies of the Killing horizon. The Killing horizon has spherical topology consisting of an inner and outer 
sphere. The event horizon is located outside of the outer Killing horizon, and the central object is a black hole. Center: 
£3 — 0.32 « e™^''~''''(a). The inner and outer Killing horizons merge in the equatorial plane, and the event horizon vanishes. 
Right: £3 — 0.33 > £™'"'~'"^(a). The Killing horizons about the origin have split into two sphere-like surfaces located above and 
below the equatorial plane, respectively. In the central and right panels, the central object is a naked singularity located at the 
Killing horizon. For all values of the parameter £3, the origin is likewise singular. 
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FIG. 5. Different shapes of the Killing horizon in the MK metric in the case £3 = — 1 for several values of the spin. Left: 
\a\ — 0.9AI. The Killing horizon has spherical topology consisting of an inner and outer spherical surface. Center: \a\ = M. 
The inner and outer Killing horizons merge at the poles. Right: \a\ — I.IM. The topology of the Killing horizon is toroidal. In 
all cases, the central object is a naked singularity located at the Killing horizon. For all values of the spin, the origin is likewise 
singular. 



A. Lorentz Violations 

When I say that a metric is of Lorentzian signature, I 
mean that 



det(5^^) < 0. 



(74) 



which is clearly negative definite everywhere outside the 
singularity (as well as the poles, = 0,7r). Similarly, 
I will study whether the metrics proposed in Section |TT] 
remain of Lorentzian signature everywhere outside their 
central objects. 



In the case of the Minkowski metric, the determinant is 
simply —1. For the Kerr metric, the determinant is 



det {g^,) ^ -J:^ sin^ 6 ., 



(75) 



B. Closed Timelike Curves 

In accordance with the no-hair theorem, the exterior 
domain of the Kerr metric, i.e., the domain outside the 
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event horizon, is causally well-behaved and free of closed 
timelike curves if |a| < M [6al- This follows from the fact 
that, for constant times t, the hypersurfaces (r, 9, </>) are 
always spacelike. This comes about because for constant 
{t,(j)), the two-dimensional metric induced on surfaces 
(r, 9) is positive definite and 



V ''(0) 



>o. 



(76) 



where C^^^ 



(0, 0, 0, 1) is the axial KiUing vector 3, HI]. 
If \a\ > M, the event horizon disappears and a naked 
singularity emerges. In this case, causality is violated 
everywhere, because any event in that spacetime can be 
connected to any other event by both a future and a past 
directed timelike curve [J,[6a| (see, also, Ref. |67|). 

Closed timelike curves may exist if the metric is no 
longer Lorentzian or if CaiC^ < 0- In the latter case, e.g. 
circles with t — const, r — const, 9 = const are closed 
timelike curves, because the vector ^^ is timelike, which 
can be thought of as the "tipping over" of light cones. In 
the Kerr spacetime with \a\ < M , closed timelike curves 
exist inside the (inner) event horizon r_ [66J. Other ex- 
amples of spacetimes with closed timelike curves include 
the van Stockum spacetime 68] , the Godel universe [63| , 
and the Gott cosmic string [70|. 

If a given metric contains a pathological region out- 
side of the central object that I identified in the previous 
section, I also calculate the largest radius at which the 
metric becomes pathological as a function of both the 
spin and the parameter C. This outermost radius serves 
as an indicator that a cutoff radius has to be introduced 
which properly excises the pathological region. 

In this section, I discuss the determinants of the QK, 
BK, and MGBK metrics to first order in the deviation 
parameters, while I discuss the full determinants of the 
MN and MK metrics. If I likewise analyze the determi- 
nants of the QK, BK, and MGBK metrics to all orders 
in the deviation parameters, I find that, if regions where 
Lorentz symmetry is violated are present, their locations 
are slightly modified, but qualitatively the same. My 
findings for the regions containing closed timelike curves 
are unchanged because of the linear form of the (0, (/)) 
elements in these metrics. 



1. QK Metric 

To linear order in the parameter Cqk, the determinant 
of the quasi-Kerr metric is given by the expression 

„3 



det(.9««) 



5ec 



2M{2M^ 



16Af2 
3Mr - Sr^ 



1 -I- 3 cos 29) 



+3r{r^ - 2Af 2) In 



r-2M 



(77) 



If Eqk 7^ 0, this determinant changes sign outside of the 
singularity at r = 2M, and hence the QK metric becomes 
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FIG. 6. Null surface and regions with Lorentz violations and 
closed timelike curves (denoted "CTCs") in the QK metric 
for a — 0.3M. Top: eqK = 0.1; Lorentz violations and closed 
timelike curves occur around the poles. Bottom: eqk = —0.1; 
Lorentz violations and closed timelike curves occur near the 
equatorial plane. 



non-Lorentzian. In addition, regions with closed timelike 
curves exist outside of the singularity. Regions of Lorentz 
violations and closed timelike curves are present around 
the poles if Eqk > and near the equatorial plane if 
fiQK < 0. In Figure IH I plot these regions for the case 
a = 0.3AI and e^^ = ±0.1. 

In Figure [71 I plot the outermost radius, at which a 
Lorentz-violation occurs, as a function of the deviation 
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FIG. 7. Outermost radius, at which a pathology occurs in 
the QK metric, as a function of the deviation parameter eqK. 
This radius depends only very weakly on the value of the 
spin. Note that I extrapolated this radius to include larger 
values of the deviation parameter; in this case, the location 
of the outermost radius is only approximate. In the range 
of the parameter —0.5 < eqk < 0.5, the outermost radius 
lies well inside of the ISCO radius (see Ref. [20|]). The dashed 
line corresponds to the naked singularity located at the radius 
r = 2M. 



parameter. At each radius, pathological regions can only 
lie on or inside of this radius. For positive values of the 
parameter eqK, this outermost radius of the Lorentz- 
violating region is located at the poles, while for neg- 
ative values of the parameter eqk, the outermost radius 
is located in the equatorial plane. This radius depends 
only very weakly on the value of the spin. In Figure [7l 
I extrapolated the location of the outermost radius to 
include larger values of the parameter eqk; this loca- 
tion is only approximate. In the range of the parameter 
—0.5 < EQK < 0.5, the outermost radius lies well inside 
of the ISCO radius (see Ref. 0). 



2. BK Metric 




X (M) 



FIG. 8. Null surface and regions with Lorentz violations and 
closed timelike curves in the BK metric for a — 0.3AI and 
B2 = ±0.1. The null surface is closed if B2 = 0.1 and open if 
B2 = —0.1. In both cases, Lorentz violations or closed time- 
like curves occur outside of the null surface and the outermost 
singularity. 



In the case of the BK metric, to linear order in the pa- 
rameter Cbk, the determinant is given by the expression 



det(gP)=-E^sin^0 [1-4(^-1 
which becomes non-Lorentzian when 

V'l - 71 > I • 



71)] 



(78) 



(79) 



For this spacetime, Lorentz violations and closed time- 
like curves occur outside of the null surface and the out- 



ermost singularity ri_+ around the equatorial plane if 
B2 y^ 0- Additional regions of Lorentz violation lie in- 
side the null surface. In Figure [H I plot these regions for 
a^O.SM andSs = ±0.1. 

In Figure [9l I plot the outermost equatorial radius, at 
which the BK metric becomes pathological, as a function 
of the deviation parameter B2 for several values of the 
spin. For positive values of the parameter B2, this ra- 
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FIG. 9. Outermost radius, at which a pathology occurs in 
the BK metric, as a function of the deviation parameter Bi 
for several values of the spin a. Note that I extrapolated this 
radius to include larger negative values of the deviation pa- 
rameter; in this case, the location of the outermost radius is 
only approximate. The dashed line corresponds to the equa- 
torial singularity located at ri,+ = 2M, which is independent 
of the spin. 



dius is equal to the radius of the singularity ri__|_ given 
by Eq. (J49]). For negative values of the parameter i?2, 
this radius is equal to the radius, at which a closed time- 
like curve is located. For larger negative values of the 
parameter Bi, I extrapolated the location of the outer- 
most radius; this location is only approximate. At each 
radius, pathological regions can only lie on or inside of 
this radius. 
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3. MN Metric 

The MN metric is of a much more complex form, and 
I will not write its determinant here explicitly. Gair et 
al. 14] , Brink Tj] , and Bambi & Lukes-Gerakopoulos [48] 
analyzed this metric in cylindrical coordinates without 
the rescaling of Section [TTl They showed the existence of 
regions with closed timelike curves around the singular- 
ity. In the form used in this paper, if ^mn t^ 0, closed 
timelike curves occur around the equatorial plane out- 
side and inside of the surface that numerically coincides 
with the Kerr event horizon. In the majority of this re- 
gion, the MN metric is imaginary. In Fig. [10] I plot these 
regions for \x\ = 0.3 and ^mn = ±0.1. 

In Figure [11] I plot the outermost equatorial radius, 
at which a closed timelike curce occurs in the equatorial 
plane of the MN metric, as a function of the deviation 



FIG. 10. Singularity and regions with Lorentz violations and 
closed timelike curves in the MN metric for x = 0.3 and Qmn = 
±0.1. Lorentz violations lie inside and closed timelike curves 
occur both inside and outside of the (singular) boundary that 
coincides to the event horizon of a Kerr black hole with the 
same value of the spin. In the black shaded region the MN 
metric is imaginary. 



parameter ^mn for several values of the spin. At each 
radius, pathological regions can only lie on or inside of 
this radius. 

If I study the determinant and the (0, (p) element of the 
MN metric to linear order in the deviation parameter, I 
find that the locations of the regions containing Lorentz 
violations or closed timelike curves are shifted slightly. 
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FIG. 11. Outermost radius, at wliicli a pathology occurs in 
the MN metric, as a function of the deviation parameter qmn 
for several values of the spin a = x^ ■ The dashed line corre- 
sponds to a Kerr black hole. 



i. MK Metric 
Evaluating the determinant of the MK metric, 

,2/1 



det [g^ 






sm 



64S2 



3a^ + 8a V + Sr-* + SegM-'r 



+4a2 (2r2 +a^)cos2e + a^ cos 46*] , (80) 



I notice that it is negative semi-dcfinitc and vanishes at 
two radii rs,±{9) if €3 < — 4r+. These radii coincide 
with the locations of the Killing horizons as found in the 
previous section, because of the relation 



detig^J')^ 



MK MK 
9tt 9d>d> 



MK\ 



(81) 



see Eq. (HT|) . Therefore, the MK metric does not contain 
any Lorentz- violating regions. 

From the g^Y" element in Eq. ([T9|) I can see that for 
given values of the radius and the spin closed timelike 
curves only occur for values of the parameter 






(82) 



rameter in Eq. ((82|) have an upper bound at 



,CTC 



r^[r^ + o^{r + 2M)] 
a2(r + 2Af) 



(83) 



corresponding to the polar angle 6 = 7r/2. Solving 
Eq. (|57)) in the equatorial plane for the parameter £3, 
I find 



^hor 



^A 



l2M3 ■ 



(84) 



It is easy to see that £3 < eg™ for all values of the 
radius and the spin. Consequently, at least in the equa- 
torial plane, closed timelike curves always lie inside the 
inner Killing horizon. Generally, I find numerically that 
the region containing closed timelike curves is always lo- 
cated inside of the outer Killing horizon. Therefore, the 
MK metric does not contain any closed timelike curves 
outside of the central object (see also Ref. [35[)- 




CTCs 
Killing Horizons 



0.5 



1 1.5 

X (M) 



FIG. 12. Inner and outer Killing horizons and regions with 
closed timelike curves in the MK metric for a — 0.3A/ and 
£3 = —0.1. The metric is regular outside of the outer Killing 
horizon, which is singular and located just outside of the null 
surface. Regions with closed timelike curves occur only inside 
the inner Killing horizon. 



No closed timelike curves exist if £3 > 0. In Figure [T^ 
I plot this region for a — 0.3M and £3 = — 0.1. I can 
see that the region containing closed timelike curves is 
located inside of the inner Killing horizon in analogy to 
the Kerr metric, where closed timelike curves lie inside 
the inner event horizon. The values of the deviation pa- 



For the determinants and (0, 0) elements of MK met- 
ric expanded to linear order in the deviation parameter, I 
find that Lorentz violations occur inside of the central ob- 
ject, while the locations of the regions containing closed 
timelike curves are shifted slightly (still remaining inside 
of the central object). 
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FIG. 13. Event horizon and regions with Lorentz violations 
and closed timelike curves in the MGBK metric for a = 0.3M 
and 7i,2 = 73,1 = 73,3 = 74,2 = 0.1. The metric is regular 
outside of the event horizon, and regions with Lorentz viola- 
tions and closed timelike curves occur only inside the event 
horizon. 



5. MGBK Metric 

In the case of the MGBK metric, the determinant to 
hnear order in the deviation parameters is, again, a very 
long expression, and I will not write it explicitly. The 
metric is regular everywhere outside of the event hori- 
zon, and regions with Lorentz violations and closed time- 
like curves occur only inside the event horizon as long 
as |7z,n| < 0.1. In Figure [131 I plot these regions for 
a = 0.3M and 71^2 — 73,1 — 73,3 — 74,2 — 0.1. Violations 
can occur outside of the horizon when |7i,„| > 0.1. 

In Table 1, I summarize the properties and pathologies 
of the various metrics. 



VI. DISCUSSION 

In this paper, I compiled several parametric frame- 
works of metrics that deviate from the Kerr solution. 
Due to the special nature of the Kerr metric in general 
relativity as a consequence of the no-hair theorem, all 
such parametric spacetimes have to differ in at least one 
of the properties of the Kerr metric, which I analyzed in 
detail. 

The MN metric as well as the QK and BK metrics 
in the appropriate limits are vacuum solutions in general 
relativity, while the MK and MGBK metrics are not. The 
MGBK metric admits four conserved quantities and is of 
approximate Petrov type D, while the other four metrics 



generally admit only three conserved quantities and are 
of Petrov type I. All five metrics are asymptotically fiat, 
which I showed explicitly for the QK and MN metrics. 
The MN metric, however, requires an appropriate coor- 
dinate transformation and a rescaling of the mass. The 
QK, BK, and MGBK metrics are designed as linear de- 
viations from the Kerr metric, while the MN and MK 
metrics are non-linear parametric deviations. 

I described the calculation of event horizons in sta- 
tionary, axisymmetric, asymptotically flat metrics using 
numerical relativity techniques. I applied this approach 
to the spacetimes I studied in this paper and showed 
that the QK and BK metrics harbor naked singularities. 
The MN metric likewise describes a naked singularity 
[34| . while the MGBK metric harbors a black hole. I 
also showed that the MK metric contains a naked singu- 
larity of spherical topology for values of the parameter 
£3 < £™^''~'=i and in the form of two disjoint spherical 
surfaces if £3 > £™^'^~''i^ gee Eq. ((66)) . I determined this 
bound analytically. If treated as a small perturbation of 
the Kerr metric, the MK metric likewise describes a black 
hole for nonzero values of the deviation parameter and all 
values of the spin. I also identified regions with Lorentz 
violations or closed timelike curves outside of the central 
objects of the QK, BK, and MN metrics, while the MK 
and MGBK metrics are free of such pathologies exterior 
to the naked singularity and event horizon, respectively. 

All of these metrics can be used for astrophysical tests 
of the no-hair theorem in either the electromagnetic or 
gravitational-wave spectrum. However, the existence of 
pathologies impacts the utility of some of these metrics 
for such tests. In order to shield the adverse effects of the 
pathological regions outside of the central object from 
distant observers, a cutoff radius needs to be defined as 
an inner boundary of the exterior spacetime, where a 
given metric is free of any unphysical behavior. Such 
a cutoff acts as an artificial horizon, which "captures" 
any matter or radiation that passes through it, i.e., any 
particle entering the domain inside the cutoff leaves the 
exterior spacetime permanently. 

Tests of the no-hair theorem with both electromagnetic 
and EMRI observations probe radii that are compara- 
ble to the location of the ISCO. In many accretion disk 
models the ISCO often marks the inner disk edge, while 
future gravitational- wave detectors will be most sensitive 
to EMRIs occuring at radii roughly in the range between 
the innermost stable orbit and 10 — 20Af . Therefore, it 
is critical for such tests that this region can be properly 
modeled. In the Kerr metric the location of the ISCO de- 
creases for increasing values of the spin and merges with 
the coordinate location of the event horizon in the limit 
a — >■ M. In parametrically deformed Kerr spacetimes, 
the location of the ISCO generally depends on the spin 
and the deviation parameters (see, e.g., Ref. [20j). Thus, 
if a cutoff radius has to be introduced in a given met- 
ric, the metric can only be used for tests of the no-hair 
theorem for values of the spin and the deviation param- 
eters for which the ISCO lies outside of the cutoff radius 
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Table 1. Properties of Kerr-like Metrics 



Metric Deviation Ricci Petrov Central 

Type Flat Type Object 



Lorcntz Closed Timelikc 
Violations'' Curves'' 



Reference 



QK 


linear'' 


yes 


I 


naked singularity 


yes 


yes 


BK 


linear'' 


no<= 


I 


naked singularity 


yes 


yes 


MN 


non-linear 


yes 


I 


naked singularity 


no 


yes 


MK 


non-linear 


no 


I 


naked singularity'' 


no 


no 


MGBK 


linear'' 


no 


D 


black hole 


no 


no 



Glampedakis & Babak (2006), Ref. [12] 
Vigeland & Hughes (2010), Ref. |16] 
Manko & Novikov (1992), Ref. [34] 
Johannsen & Psaltis (2011), Ref. [35] 
Vigeland, Yunes, fc Stein (2011), Ref. [18] 



'^ Lorentz violations and closed timelike curves refer to pathologies outside of the central objects. 

'' I analyzed the properties of the spacetimes with a linear deviation from the Kerr metric only for small values of the respective 

deviation parameters, c.f., Eq. JSj. These properties may be different for larger deviations. 
'^ Ricci-flat if a = 0. 
'' Black hole for small deviations at linear order. 

TABLE I. Properties of the Kerr-like metrics: The quasi-Kerr metric (QK) modifies the quadrupole moment of the Kerr metric 
after the quadrupole moment of the Hartle-Thorne metric that is a solution of the vacuum Einstein equations for small values of 
the spin and the quadrupolar deviation. The bumpy Kerr metric (BK) augments the Kerr metric with arbitrary small multipolar 
distortions, with a functional form chosen to insure that it remains a vacuum solution of the linearized Einstein euqations in 
the zero spin limit. The metric proposed by IVIanko & Novikov (IVIN) is a Ricci-flat generalization of the Kerr metric with 
arbitrary multipole moments; it is an exact solution of the Einstein equations. The metric proposed by Johannsen & Psaltis 
(MK) deviates from the Kerr metric in non-linear form and introduces a set of polynomial modiflcations to the spacetime. 
The modified-gravity bumpy Kerr metric (IVIGBK) is an extension of the BK metric designed to admit an approximate fourth 
constant of the motion. 



(see discussion in Refs. [20, Is^)- This implies that the 
BK and IVIN metrics (as well as the QK metric by con- 
struction) can only be used for moderately spinning black 
holes as long as the ISCO lies outside of the cutoff radius. 
For these three spacetimes, I calculated the outermost ra- 
dius, at which the spacetime is no longer well behaved, 
as a function of the spin and the deviation parameter. 
In each case, a cutoff radius has to be introduced so that 
the pathological regions are excised. In practical applica- 
tions, it is often convenient to choose a spherical cutoff, 
which is located just outside of this outermost radius. 
However, more sophisticated choices of a cutoff radius 
which depend on the polar angle 6 are also possible. 

The IVIK metric is regular everywhere outside of the 
naked singularity located at the Killing horizon for all 
values of the spin \a\ < M. In the region of the pa- 
rameter space where an ISCO exists, the ISCO is always 
located outside of this surface and, as shown in Ref. [33 , 
the ISCO only coincides with the Killing horizon in the 
Kerr case if a = M and £3 = 0. Therefore, if needed, 
a cutoff can always be chosen at a radius between the 
naked singularity and the ISCO (a similar property also 
holds for the circular photon orbit; see Ref. [35]). 

In several settings of accretion flows around black 
holes, the presence of an event horizon affects only 
marginally the characteristic of the flow itself or its obser- 
vational appearance. Inside the ISCO (i.e., in the plung- 
ing region) the plasma attains highly supersonic inward 
velocities and gets causally disconnected from the hydro- 
dynamics of the material outside the ISCO. Imposing a 
cutoff, therefore, only has a very minor effect, and the 
MK metric provides a useful framework for electromag- 
netic tests of the no-hair theorem for all values of the 
spin |a| < M [351]. In other situations, however, the pres- 



ence and properties of the event horizon modifles sig- 
niflcantly the accretion flow outside the ISCO. This is 
especially true for spinning black holes with significant 
magnetic flux threading the horizon, potentially arrest- 
ing the flow, and launching powerful jets. In such cases, 
(e.g., in fully relativistic magnetohydrodynamic simula- 
tions), more care must be taken to ensure that flelds and 
material in the simulation do not come into contact with 
a singularity. I will analyze the utility of the MK metric 
in this scenario in a future paper. 

The MGBK metric is regular everywhere outside of 
the event horizon as long as it remains a small perturba- 
tion of the Kerr metric. In this regime, this metric can 
likewise be used for electromagnetic tests of the no-hair 
theorem. However, since this metric is of Petrov type D, 
it admits an approximate fourth constant of the motion, 
a Carter constant. This property is useful to construct 
model waveforms, and the MGBK metric is, therefore, a 
useful framework for gravitational- wave tests of the no- 
hair theorem [l^ [l^ . 
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Appendix A: Coordinate Transformations in the 

MN Spacetime 

I begin with the line element as given by Ref. [3J] (see, 
also, 14]) 

dx"^ dy'^ 



1 1 



(Al) 



where 



J MN — c , 



= 2fce-"*^ 



C 



4fcQ! 



g2r ^ g2r' 



A 1 - a2 ' 



(.T2-l)(l-a2)2' 

A = {x^- 1)(1 + Mi;)2 - (1 - y^){v - uf, 
B = [x + l + {x- 1)™]^ + [(1 + y)u + (1 - y)vf, 
C = {x^ - 1)(1 + w)[w - u - y(u + w)] 
+ (1 - j/^)(w - u)[l + uw + x{\ - uv)\, 



1 



x^- 1 



2 x^ — y^ 



^Pi Pi) 



J=0 



X- j/+(-l)^~'(a; + j/) 



Pi -2 



u = —a exp 



-^''(-i+E^^ 



V = a exp 



2/3 1+5: 



/=0 



{-lY-'{x + v)Pi 



i?= v^a:^ + y2 _ 1^ 






)' ^"(^)-2i!(i:)"(^^-^)^ 



(A2) 

In these equations, fc, a, and /3 are free parameters, which 
determine the set of multipole moments of the spacetime. 
Note that there is a minor discrepancy in the functions 
u and V (a and h in their notation) between the original 
form by Ref. 34] and the one used by Ref. [1J|. The 
difference, however, is small, and both versions satisfy the 
vacuum Einstein equations. I will use the version quoted 
by Gair et al. [14[ in this paper. Following Ref. [ij], I 
define 



-M + ^/^P 



k = M 
P = q 



a 

1 + ^2 



= ^/AP - a? 



(A3) 



where M and a are the mass and the spin, respectively. 
This metric can be mapped to cylindrical coordinates: 



p = k\J x'^ — ly 1 



y 



z = kxy 



(A4) 



X 



V 



with inverse 

1 - 

2fc . 
1 - 

2fc . 
The MN metric then becomes 



v/p2 + (z + fc)2 + ^p2 + (z - fc)2 
v/p2 + (z + fc)2 - ^p2 + (z - fc)2 



(A5) 



c^sL = -/MN(di-c^d<^)'+/MN [e''^(V + rfz') + P'd^^\ . 

(A6) 
In spherical coordinates, I define [3J] 



r = fca; + M, cos = y, 

and the MN metric is as given in Eq. (fT6|) . 



(A7) 



Appendix B: Other Forms of the MGBK Metric 

Certain choices of the free parameters of the MGBK 
metric allow for a mapping to other alternative theories 
of gravity. See Ref. [18j for a detailed discussion. Here I 
calculate the event horizons of the black hole metrics in 
Einstein-Dilaton-Gauss-Bonnct and Chern-Simons grav- 
ity using the condition in Eq. (j47p . 



1. Linearized Kerr 

First, however, I consider the Kerr metric, linearized in 
a^, which can be regarded as a "deformed Schwarzschild" 
spacetime. Expanding the metric element 



A 



(Bl) 



to C(a^) and defining h^ by subtracting the correspond- 
ing element of the Schwarzschild spacetime, 



Schw 



!s{rr 



where 



I find 



fsir) ^ 1 



2M 






Eq. dill) becomes 



[cos^ e ~ fsir)-^] 



(B2) 



(B3) 



(B4) 



fsirYh: 



2udS 



:[fs{r) cos^e- I] ^0. (B5) 
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When solving this equation, I need to bear in mind that 
my analysis is done to 0{a^), and so the horizon will be 
shifted by this amount [recall Eq. (|44p ]. I insert 



(B6) 



into Eq. (jB5[) . linearize in a^, and find 
2Ma2(l + 4A) = 0. 



(B7) 



Thus, A — —1/4, and the horizon of my linearized Kerr 
spacetime is at 



rH = 2M 1 



1 a' 



2M 



la' 



(B8) 



4M2y '■' 2Af' 

This is of course is just the exact solution for the Kerr 
horizon, r+ = AI + y/AP — a^ , expanded to 0{a'). 

2. Einstein-Dilaton-Gauss-Bonnet Gravity 

Yunes & Stein [40| looked at black holes in a class of 
gravity theories described by Lagrangians modified from 
the standard Einstein-Hilbert form by scalar fields cou- 
pled to quadratic curvature invariants. They find a set 
of such theories that admit black hole solutions, though 
deformed from the usual Kerr form. Focusing on non- 
rotating solutions, they find that for a class of such solu- 
tions the relevant component of the metric deformation 
is given by 



EDGE 



"3 



M 52 IVP 

1 + — + ^ 
r or 



KAPr'fsir)' 
2Af3 



'4^-'-^^)m 



The constant a^ is the theory's coupling constant, and 
K = l/(167r). Note that this spacetime can also be re- 
garded as a "deformed Schwarzschild" case, with a form 
for 7i(r) that is somewhat different than that used in 
Ref. [ii]. 

I use this /i™°^ in Eq. gT]), and use g^chw ^ j^,,^j,y 
My solution should now take the form 



rn =2M{1 + Xa3). 



(BIO) 



Inserting this expression into Eq. (j47l) and linearizing in 
as, I find 



A-^^^=0, 



80 M^K 



and thus 



rn = 2M [1- 



49 as 

SOAf^K 



(Bll) 



(B12) 



This expression agrees with the text following Eq. (12) 
of Ref. [4fl]. 



3. Chern-Simons Gravity 

As a final example, I consider slowly rotating black 
holes in dynamical Chern-Simons gravity, a case analyzed 
by Yagi, Yunes, & Tanaka 6J|. The relevant metric com- 
ponent in this case is given by 



^cs 



'fs{r) 



[cosH - fs\r)] + Cx' 



M^ 



r'Psir) 



201 
896 



fs{r)ni{r)P2{cos0) 



25 M 
384T 



n2{r) 



(B13) 



Here, Cx^ are a combination of coupling parameters that 
determine the strength of the Chern-Simons modification 
to gravity to this order, P2 = (3 cos^ 9—l)/2 is a Legendre 
polynomial, and ni.2(r) are polynomials in (M/r), which 
are given in Ref. [40]. 

In this case, I expect the horizon to be at 



(B14) 



r«=2Af (1 + Ai— + A2Cx^ 



Building the horizon equation (|T7|) . inserting this form 
for rn, and linearizing in a' and Cx^i I find 



a^ /I 
Af2 I 4 



+ '^1+^^ (5™ 



+ A2 = 0. 



(B15) 



This means 



1 a2 „ 915 

^--2*^1 1-4M^-^^ 57344 



2Af 



915 



<X^M. 



(B16) 



28672 
This agrees with Eq. (55) of Ref. |64,], linearized in a' 



Appendix C: Numerical Methods for Finding Event 
Horizons 



Here I describe two numerical methods for solving 
Eq. P5)) in order to construct the event horizon. 



1. Spectral Method 

The following technique only works if the horizon's ge- 
ometry can be described as a "Strahlkorper" [53] ■ A 
Strahlkorper is a figure that encloses the origin and is 
only intersected once by any ray from the origin. This is 
the case for mild deformations from Kerr, but may not 
be the case if the deformation is more severe. A horizon 
of unusual topology (e.g., disjoint horizons on the z-axis, 
or a torus) is certainly not a Strahlkorper, at least not 
about "r = 0". It may be possible to adapt this tech- 
nique to such situations by changing the origin (e.g., to 
the middle of a disjoint horizon, or to the center ring of 
a torus). 
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Bearing in mind my boundary conditions [H{0) = Hk, 
H{tt/2) = i?eq, dH/d9 = at 6* = 0,7r/2], I proceed by 
writing 



N 



H{9)=Y,'^nPnicOSt 



(CI) 



ra=0 



where P„ are the Legendre polynomials and where iV is a 
suitably chosen truncation of the sum over them. Solving 
for H{0) then means solving for the expansion coefficients 
a„. 

All of the spacetimes that I consider are reflection sym- 
metric about 9 = 7r/2, so I need only include even values 
of n. I then modify the expansion to 



N 
H{9)=^a2nP2n{cOS( 
n=0 



(C2) 



I choose A^ — 1 angles evenly spaced between 9 = and 
9 = it/2. By enforcing Eq. (H51) at these A^ — 1 angles, 
plus the solutions at 9 = and 9 = 7r/2 from Eq. ([57)) . I 
have a total of A" + 1 equations for the A^ + 1 coefficients, 
completely specifying the horizon function H{9). 



2. Finite Difference Method 

A crude but robust method for solving the horizon 
equation is to use finite difference to construct the deriva- 
tive and solve the equation. Cover the angular sector 
with A^ segments of width 59, so that there are a total of 
A' + 1 angles from 6*0 = to ^at = 7r/2. Let Hi denote 
the solution at angle 9i. Then, the horizon equation (j43|) 
becomes 



= g'-^{H,,9,)+g'''{H,,9,) 



Hi-i — Hi 



(C3) 



The solution is known at i = (0 = 0); from there, step 
to i = N {9 = 7r/2). Empirically, I have found that this 
method agrees with the spectral method extremely well if 
there is a null surface that encloses the origin and passes 
through 6* = and 6* = tt/2. 

It may be possible to generalize this method to handle 
different null surface topologies. Since I have not iden- 
tified any parameter choices that unambiguously yield 
such unusual null surfaces, I have not tested this. Note 
that the null surfaces in the MK metric, where they exist, 
always have spherical topology, while the Killing horizon 
can have spherical, disjoint, or toroidal topology. 
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